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Abstract 

The analytic calculation of the three-loop QCD /3-function and anoma- 
lous dimensions within the minimal subtraction scheme in an arbitrary 
covariant guage is presented. The result for the /3-function coincides with 
the previous calculation pi. 
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The calculation of the one-loop /3-function in Quantum Cromodynamics 
(QCD) has lead to the discovery of the asymptotic freedom in this model and 
to the establishing QCD as the theory of strong interactions |^]. The two-loop 
QCD /3-function was calculated in ||. The first calculation of the three- loop 
/3-function in QCD was done in Q . Since then the perturbative calculations of 
physical quantities in QCD reached the (next-next-to-leading) order of per- 
turbative expansion in the strong coupling constant a s for the total cross section 
of e + e~~ annihilation into hadrons Q and for the deep inelastic lepton-hadron 
scattering || . The phenomenological application of these calculations involves 
the three-loop approximation of the /3-function. The three-loop /3-function is 
also necessary for the investigation of the scheme-dependence problem in QCD 
and infrared behavior of the coupling constant That is why it is im- 

portant to have an independent check of such a complicated calculation as the 
calculation of the three-loop QCD /3-function ]|]. 

In this paper we present original three-loop results for QCD anomalous 
dimensions in an arbitrary covariant gauge and confirm the three-loop result 
for the /3-function of the work J1J]. 

Throughout the calculations we use the dimensional regularization |p and 
the minimal subtraction (MS) scheme |J. The definition of the /3-function is: 

/3(a) = ^ = -5>a*+ 2 , (1) 

here a = g 2 /lQir 2 = a s /4ir and g is the strong coupling constant in the stan- 
dard QCD lagrangian. To calculate the /3-function we need to calculate those 
renormalization constants of the lagrangian which determine the renormaliza- 
tion constant of the charge. It was convenient for us (for technical reasons) to 
choose the following three renormalization constants: Z\ of the quark-quark- 
gluon vertex and correspondingly Zi of the inverted quark propagator and Z3 
of the inverted gluon propagator. Note that the authors of the work [jlj] used 
another combination: ghost-ghost-gluon vertex and correspondingly inverted 
ghost and gluon propagators. So our calculation provides a really independent 
check of the three-loop /3-function. 

Renormalization constants within MS-scheme do not depend on dimensional 
parameters (masses, momenta) ]l0| and have the following structure: 



00 



Z(a) = l + V^-p, (2) 



n=l 



where e is the parameter of the dimensional regularization, the dimension of 
the space-time being D = 4 — 2e. The renormalization constants define corre- 
sponding anomalous dimensions: 

2 d\ogZ(a) dzW(a) 

7(a) = " = - a -&r- (3) 

It is convenient to express the /3-function through the corresponding anomalous 
dimensions: 

/3(a) = a [ 73(a) + 2 72 (a) - 2 7l (a) ]. (4) 
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The calculation of the renormalization constants within MS-scheme can be re- 
duced to the calculation only of massless propagator diagrams by means of the 
method of infrared rearrangement [11]. In the case of the quark-quark-gluon 



vertex it means in particular that we can safely nullify the gluon momentum 
reducing the calculation of the Z\ to the propagator diagrams. 

One way to find the constants Zi is to obtain them as the sums of countert- 
erms of individual diagrams. This implies application of R-operation to each 
individual diagram which is an enormous amount of work in the three-loop 
case. This way was used in [0. But now we have en efficient program Mincer 
12] written for the symbolic manipulation system Form 12]. This program 



computes analytically three-loop massless propagator diagrams within dimen- 
sional regularization. So we are able to obtain renormalization constants at the 
three-loop level straightforwardly from the multiplicative renormalizability of 
Green functions: 

r Renormalized(a) = Z(d) TBare(aB), (5) 

where as is the bare charge. 

Thus we compute with the program Mincer the unrenormalized three-loop 
one-particle-irreducible quark-quark-gluon vertex and inverted quark and gluon 
propagators. Having the two-loop bare charge we determine the necessary three- 
loop constants Z\, Z2 and Z^ from the requirement that the poles in e cancel in 
the r.h.s. of eq.(|5j). The computations were done in the covariant gauge with an 
arbitrary gauge parameter £ to check directly the gauge independence of the (3- 
function. So the diagrams were run with the gluon propagator (g^ — £^p^) / Q 2 ■ 
(Note that the computations in were done in the Feynman gauge £ = 0.) 

The results of our three-loop computation of the anomalous dimensions in 
the MS-scheme are: 

7l (o) = +a[C F + C A + i{-C F -l/AC A )] 

+« 2 [ - \cl + y CfCa + %p A " Cpnf " T2° Anf 
+ t{-\c F c A - y q c 2 a) + e(\c F c A + \c\) ] 

+a 3 [ \C% + ClC U(12C 3 ~ l -f) + C F Cl(-^C 3 + ^f) + Cl(\cs + ^) 
+ \c\n f + C F C A n f (6C3 - ^) + Cln f (- 9 - C 3 - ^) + \c F n) - ^C A n) 
+ Z(C F C A (- S -ts - ^) + C\(- A C3 - ^) + ^C F C A n f + \c\n f ) 

+ e(c F cl { \c, + § ) + ci(| C3 + 1 )} + e(-^c F cl - LcD ]. (6) 

72 (a) = +a [ C F + a~C F ) ] 

+a 2 \-\ C F + y CfCa -° Fn f + t(-\c F C A ) + e{\c F C A ) ] 
+a 3 [ \CI + a^(12Cs - ^) + C F C\{- 1 ^ + ^) + 3 -C F n f 
- ^-C F C A n f + \c F n) +^C F C 2 A (~C3 ~ ^) + y^C^) 
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+ e(C F C\( 3 -C3 + §)) + e{-y Q C F C\) ], (7) 
73(a) = +a [ --C A + -n f + Z{--C A ) } 

+a 2 [ ~ f C\ + 2C F n f + \c A n f + ^C A ) + e{\c\) ] 

o 4051 5 875 

+« 3 PK-Cs - W ) - C^n/ + CjpCxn/ClJ^ + -) + C\n } (-%, + —) 
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+ C(CI(-|C3 ~) + C\n f ) + e 2 (Ci( Afc + g)) + e 3 (-|ci) ]. (8) 



Here Cp = {v? c — l)/2n c and C A = n c are the Casimir operators of the 
defining and adjoint representations of the color group SU (n c ), n/ is the number 
of active quark flavors, £3 is the Riemann zeta-function (£3 = 1.202056903 . . .), 
£ is the gauge parameter. 

Substituting the obtained 7- functions to the eq.(Q) we find the desired (3- 
function: 

/3(a) = - a \+^CA-\n f \ 

-a s [^C A -2C F n f -™C A n f ] 

4r 2857^ 3 ^ 2 205^ „ 1415^ 2 11 2 79^ ^, 

"° hgT^ + C ^ " l8 CFCAnf ~ ~54T CATlf + ~~g f + m Ca ^ } 

This result coincides with the previous calculation ]l]]. 

The /3-function is gauge independent as it should be within the MS-scheme 
and does not contain £3. The cancellation of the £-terms is a strong check of 
the validity of the calculations. 

We repeat the result once more for the particular case of QCD (n c = 3): 

,,, . a, 2 . o. 38 . ji ,2857 5033 325 o x . . 
/3(a) = -a 2 (ll - -n f ) - a 3 (102 - -n f ) - a 4 (— —n f + — n}). (10) 

It is interesting to note that the poles of individual three-loop diagrams of 
the quark-quark-gluon vertex and the the total renormalized three-loop quark- 
quark-gluon vertex (which we have also calculated but do not present here 
because of the cumbersome result) contain the quartic Casimir operator (see 
e.g. (l4|])- But it cancels in the anomalous dimension 71 and correspondingly in 
the /3-function which are expressed through the quadratic Casimir operators C F 
and C A only. Since the quartic Casimir operator is present in the finite part of 
the three-loop quark-quark-gluon vertex it can enter the four-loop /3-function. 



But although existing techniques 11], [12| are sufficient to compute the four 



loop /3-function analytically, in practice it is a very difficult task because of 
enormous amount of diagrams to be calculated. 
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